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In quantum wires, the hyperfine coupling between conduction electrons and nuclear spins can lead 
to a (partial) ordering of both of them at low temperatures. By an interaction-enhanced mechanism, 
the nuclear spin order, caused by RKKY exchange, acts back onto the electrons and gaps out part of 
their spectrum. In wires with two subbands characterized by distinct Fermi momenta kpi and 
the nuclear spins form a superposition of two helices with pitches ir/kFi and n/kF2, thus exhibiting 
a beating pattern. This order results in a reduction of the electronic conductance in two steps upon 
lowering the temperature. 

PACS numbers: 71.10.Pm, 75.30.-m, 73.22.-f, 75.75.-c, 31.30.Gs 



I. INTRODUCTION 

Being a hallmark of topological states of matter as well 
as a versatile platform for engineering quantum com- 
putation devices, helical and quasi-helical electron sys- 
tems have attracted much interest in recent years. Be- 
sides their prominent formation as edge states in Quan- 
tum Spin Hall samples^ (quasi-) helical Luttinger liq- 
uids can for instance also be generated by applying a 
magnetic field to a Rashba spin-orbit coupled quantum 
wire^ When brought into contact with an ordinary su- 
perconductor, a helical Luttinger liquid can turn into a 
topological superconductor with zero energy Majorana 
bound states at the two ends of the wire^— Helical Lut- 
tinger liquids are however also interesting for their own 
sake and have for instance applications as spin filters. 2 

In addition to spin-orbit based proposals, it has been 
realized that quasi-helical Luttinger liquids may also 
emerge due to helical magnetic fields. The helical field is 
thereby fully equivalent to the combination of a homoge- 
neous magnetic field and Rashba spin-orbit interactions. 6 
The interplay of such helical fields (and more generally 
oscillating fields, for instance generated by nanomagncts) 
with the usual spin-orbit coupling furthermore allows to 
extend the topological phase diagram of Rashba spin- 
orbit coupled quantum wires to exotic phases beyond 
Majorana bound statesi 7 -^ 

As an alternative to an applied helical magnetic 
field, the intrinsic hyperfine coupling between electrons 
and nuclear spins can lead to the spontaneous forma- 
tion of a quasi-helical Luttinger liquid at low tempera- 
tures. As has been argued in the case of a single sub- 
band quantum wire, the hyperfine coupling results in 
a Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction 
that diverges at momentum 2k p due to low energy elec- 
tron backscattering. 9 ' 10 This in turn induces a helical or- 
der of the nuclear spins, which feeds back to the electrons 
as a magnetic Overhauser field. The latter gaps out half 
of the electronic degrees of freedom, thereby turning the 
quantum wire into a quasi-helical one. Importantly, the 
electronic gap, strongly enhanced by electron-electron in- 
teractions, opens exactly around the Fermi level with- 
out any fine tuningi&i - We expect the helical (nuclear) 



spin order to be observable not only in hyperfine cou- 
pled quantum wires, but more generally in systems that 
fall into the class of weakly-coupled, electronically one- 
dimensional Kondo lattice systems of finite size in the 
so-called RKKY liquid phase when both the Kondo en- 
ergy scale k,B Tk and the direct exchange between the 
Kondo lattice spins are smaller than the RKKY interac- 
tion (the nuclear spins are one possible example of such 
a Kondo lattice) *ii The helical order of the nuclear spins 
in quantum wires is furthermore intimately related to a 
proposed ferromagnetic order of the nuclear spins in two- 
dimensional electron gases 

The quasi-helical state resulting from gapping out 
parts of the electronic spectrum by either a combination 
of a homogeneous magnetic field and spin-orbit interac- 
tions, a helical magnetic field, or a spontaneous ordering 
of the nuclear spins, can in some ways be understood 
as a helical Luttinger liquid (the spin of the remaining 
gapless electronic modes is approximately locked to their 
direction of motion). Their physics is however even richer 
than the one of an ideal helical liquid. This has been dis- 
cussed in terms of the spectral properties of the conduc- 
tion electrons^ For a better distinction from ideal heli- 
cal Luttinger liquids, the quasi-helical systems resulting 
from gapping out parts of the electronic spectrum have 
been dubbed spiral Luttinger liquids. Helical Luttinger 
liquids can in fact be understood as a special subclass of 
spiral Luttinger liquids with charge and spin Luttinger 
parameters K c = 1/K S J£ The study of spiral Luttinger 
liquids is consequently of specific interest as new phenom- 
ena (beyond helical Luttinger liquid physics) can arise. 
Since furthermore a considerable number of experiments 
and theoretical proposals are based on spiral Luttinger 
liquids mimicking helical Luttinger liquids, the under- 
standing of spiral Luttinger liquids is also of practical 
relevance. 

In this work, we discuss how stable the spontaneous 
formation of helical nuclear spin order in a quantum 
wire is to the presence of multiple subbands. Indeed, 
multi-subband quantum wires are characterized by dif- 
ferent Fermi momenta kp% for the different bands. One 
may thus expect low-energy electron backscattering at 
any combination kpi + kpj, and it is not obvious how the 
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nuclear spins order in this case. To analyze their inter- 
play with the electrons, we specifically focus on a quan- 
tum wire with two subbands that may either correspond 
to the lowest two bands of a single wire, or alternatively 
emerge from two sufficiently coupled parallel quantum 
wires with a single occupied band each. In the latter case, 
the two subbands are the symmetric and antisymmetric 
orbitals shared between the two wires. After deriving the 
effective one-dimensional model in Sec. [Ill we list an d an- 
alyze the possible orders of the nuclear spins in Sec. IIIII 
We then discuss that the nuclear spins form two super- 
imposed helices at momenta 2kpi and 2kp2, but not at 
kpi + kp2, and analyze the resulting low energy theory 
in Sec . 11 VI in a self-consistent fashion. In Sees. W\ and IVI[ 
we finally turn to the onset and stability of this order. 



II. THE MODEL 

For concreteness, we specialize to a model of a single 
quantum wire defined in a two-dimensional GaAs elec- 
tron gas by virtue of electrostatic gates. The chemical 
potential, tunable by an electrostatic gate, is chosen such 
that only the lowest two subbands of the wire are (par- 
tially) filled. Nevertheless, our findings remain qualita- 
tively valid for Carbon nanotube samples or the above 
mentioned parallel wires. 

As a central ingredient for our discussion, the electrons 
interact via screened Coulomb interaction. The screening 
may for instance be due to mirror charges in the gate elec- 
trodes of the sample and is characterized by a screening 
length larger than the typical length scale associated with 
the transversal confinement. We neglect the weak spin- 
orbit interaction (note that the spin orbit length is much 
larger than typical Fermi wavelengths in GaAs wires and 
Carbon nanotubes) i 14 i 15 Expressing the electron Hamil- 
tonian in terms of right-moving and left-moving particles 
in the two subbands, the Coulomb interaction gives rise 
to various matrix elements for electrons close to the four 
Fermi points. If there were no nuclear spins, the electrons 
could be described as two spinful and gapless Luttinger 
liquids whose velocities and Luttinger liquid parameters 
are renormalized by a number of density-density type 
interactionSi 16 ' 17 For simplicity, we will, however, only 
keep track of charge density interactions since the lat- 
ter are typically larger than all other electron-electron 
couplings, resulting in strongly renormalized charge ve- 
locities and fairly unrenormalized spin velocities j^ r 21 i 23 
We will come back to the small effect of spin density 
interactions on a spin ordered state in Sec. |Vl 

The nuclear spins in the wire are modeled by a Kondo 
lattice type Hamiltonian. They interact with each other 
via dipolar and quadrupolar interactions. Most impor- 
tantly, they are also subject to a hyperfine coupling to 
the electrons. As will be discussed below, the hyper- 
fine coupling gives rise to an RKKY interaction between 
the nuclear spins** 22 Similar to single band wires---- this 
RKKY interaction overrules any direct coupling between 



the nuclear spins by orders of magnitude, and we will 
thus neglect the latter in the remainder. 

The Hamiltonian is derived by first linearizing the 
spectrum around the Fermi points. Choosing the wire 
to be aligned along the z-axis, we decompose the anni- 
hilation operator for an electron of spin a =J, i in band 
j = 1,2 as c ja {z) = e ikF i z R ja (z) + e~ ikF J z L ja (z), such 
that e ikF ' z R jt(T (e~ lk "^ z L jt<7 ) annihilates a right-moving 
(left-moving) particle {kpj denotes the Fermi momentum 
in band j). This yields 

H = H C + H cn , (la) 

He = Ysf dzV ^ (RU- id *) R i*- L U- i9 *) L J*) 

U f 

+ — / dz p to t Ptot , (lb) 
H en = J2A S i -I i , (lc) 

i 

where vpj is the Fermi velocity of band j and U de- 
notes the local interaction for the charge density ptot = 
V , . T | lf iz R ja + L\ a L ja ). The hyperfine coupling be- 
tween the electron spin Si and the nuclear spin Jj at 
site i is Aq, where i runs over all sites of the nuclear 
spin lattice within the support of the electronic wave 
functions. This includes typically N± ~ 50 x 50 sites 
in the transversal directions (the nuclear spin lattice is 
thus three-dimensional, see Fig. [TJ) We finally assume 
that the wire has a finite length L, and that the lattice 
constant of the nuclear spins is a. 

Since the Fermi energies Epj in the two bands are 
typically much larger than the intrinsic energy scales of 
the nuclear spins as well as the coupling between nu- 
clear spins and electrons, Epj ^> Aoj 10 ' 25 the dynamics 
of electrons and nuclear spins decouple. The former me- 
diate an effectively instantaneous RKKY interaction for 
the latter, while the nuclear spins act as a static magnetic 
Overhauser field for the electrons (in the disordered state, 
this Overhauser field vanishes). The separation of scales 
largely simplifies the following analysis. 

A. Effective one-dimensional model 

Because the electrons behave quasi one-dimensional, 
the coupling between nuclear spins and electrons is 
roughly the same for all spins within a given transver- 
sal section of the quantum wire. We can thus derive 
an effective quasi one-dimensional hyperfine Hamiltonian 
following Ref. Idj. To this end, we note that the wave 
function of an electron can be written as a product of the 
longitudinal and the transversal state. For a particle at 
site i — i±), this implies = <8)\i±). The more 
appropriate basis for the electrons, however, is defined 
by the two subbands. Denoting the two transversal wave 
functions corresponding to the two occupied subbands 
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FIG. 1: Panel (a) depicts a quantum wire with two subbands 
defined in a two-dimensional electron gas by a harmonic con- 
fining potential. The subbands have transversal wave func- 
tions ^j{x,y) — ipj(x)<giipi(y), and therefore exhibit different 
^-inversion symmetry with respect to the wire axis. A sketch 
of the nuclear spin lattice in the wire is given on the left- 
hand side of panel (b) (for visibility reasons, not all spins are 
shown). Within a given cross-section, the nuclear spins align 
ferromagnetically, while they form helical orders along the 
wire direction (see Sec. IIV|) . The right-hand side of panel (b) 
shows the effective description of the nuclear spins as helical 
magnetic fields. 



j = 1,2 as tyj(x,y), we find that the annihilation oper- 
ator for an electron of spin a at site i can be expressed 
as 



ia ~ E ' 



(2a) 



N ± / dx dy (x, y) $* ± (x, y) , (2b) 



C 



where $>ij_(x,y) is the Wannier orbital associated with 
i± (chosen such that Cji ± is a real number). Within the 
support of the electronic wave functions, Cji ± is of order 
one, while it vanishes outside the support. We can now 
express the hyperfine coupling in the two subband basis 
as 



_ Aq ^ q j 

n ~ Nj_ JJ'' l n ij'>« 



E 

a0 



Ijf,i\\ ~ E^ Cjij x Ij 



(3a) 
(3b) 
(3c) 



where <x denotes the vector of Pauli matrices. 

The terms j = j' describe how the spin of electrons in 
band j couples to the nuclear spins within a given section 
in of the wire. Since (Cji ± ) 2 ~ 1 across the section, the 
electron spin in fact rather couples to one big nuclear 
spin 



"'"Mil ~~ /XClij_) 2 Ii n i± ~ y^(C2^) 2 Jt 



(4) 



than to N± individual nuclear spins. The interaction of 
an electron with a localized spin also allows for inter- 
band scattering processes. Since the two bands either 
correspond to the lowest two subbands of a single quan- 
tum wire or the symmetric and antisymmetric orbitals 
shared between two wires, the transversal wave functions 
\l/i(x, y) and ^2(x,y) have different inversion symmetry 
with respect to the wire axis, see Fig. [T] We choose the 
first orbital to be the symmetric one, 



$i(x,y) =+V 1 (-x,y) , V 2 (x,y) = -y 2 (-x,y) ■ (5) 

The different inversion symmetry of the initial and fi- 
nal orbital of an interband spin scattering event implies 
that these processes couple to an inversion antisymmetric 
combination of the nuclear spins, 



= E \Cu ± C 2i± \sgn{i± x ) Ii^ 



(6a) 
(6b) 



With these definitions, the effective hyperfine Hamilto- 
nian reads 



A 



11,111 T J22,i 



' J S,i|| 



(7) 



E ( Sl ^H + S 21,i||) ' /a 



The electrons thus couple to either a symmetric or anti- 
symmetric combination of all nuclear spins within a given 
section of the wire, depending on whether the spin scat- 
tering event changes the transversal wave function or not. 
These effective spins have the size \I s ,i« \, l-^Mii I ~ IN±, 
where / is the length of the individual nuclear spins. Be- 
cause N± is typically of the order of (a few) thousand, 
these effective nuclear spins can be treated semiclassi- 
cally, such that for instance Kondo correlations can be 
neglected. The reduced coupling strength Aq/N± finally 
reflects the spread of the electronic wave function across 
the wire section. 
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fc, w, J, 

FIG. 2: Fermionic diagrams contributing to the (x and y com- 
ponents of the) spin susceptibilities Xij- Solid lines denote a 
particle in band 1 at, e.g., momentum k + q, frequency uj and 
spin up, dashed lines denote a particle in band 2 (note that 
the Luttinger liquid approach yields the density-density inter- 
acting versions of these diagrams). The dotted lines denote 
the momentum q emitted/absorbed by the nuclear spins. 



III. HELICAL ORDER OF THE NUCLEAR 
SPINS: POSSIBLE SCENARIOS 



Like for the single subband quantum wire^S we use 
the separation of scales between the nuclear spins and 
the electrons to step by step derive a self-consistent de- 
scription of the coupling between electrons and nuclear 
spins. We start from two gapless Luttinger liquids and 
a disordered bath of slow nuclear spins. By virtue of 
the hyperfme coupling, the electrons mediate an RKKY 
interaction for both nuclear spin superpositions I s and 
J a . As detailed in Appendix lAl the RKKY interaction is 
described by the Hamiltonian 




J, 



a,q 
> k 



k F i + k F2 
+■ 




**q 



FIG. 3: RKKY exchange interactions J a and J a for <?22 < 
<?12 < 311 and at finite temperature. 



T a P\ 

\q^±2k F1 

jal3\ 

s,Q \qa± 2 k F2 
U ^l\q^±(k F1 +k F2 ) 



<W \qT2k F1 \ 2911 - 2 , 
<W |<?T2fc F i| 2922 - 2 , 

~>S a0 \qT(k F i + k F2 )\ 2912 - 2 



(9a) 
(9b) 

(9c) 



We note that the exchange couplings are diagonal in spin 
space because the Hamiltonian preserves the total spin. 
At finite temperatures, the divergences turn into sharp 
dips, whose depth is controlled by the temperature T and 
the exponents <7-y , and whose width is of the order of the 
thermal momentum, see Fig. |3] 

In order to gain energy, the slow nuclear spins will 
now orient into the minima of the susceptibilities. As 
discussed in Refs. and [l(|, this results in helically or- 
dered states at the momenta corresponding to the dips 
of the susceptibilities. In a two-subband quantum wire, 
the presence of three different combinations of Fermi mo- 
menta k F i + k F j could a priori result in a superposition 
of the three helices, 



1 / J a * J a ' \ 

#rkky = (/"-<? jvf T L + ^-q jfi T L) > 



(8) 



where the exchange interactions are determined by the 
static part of the respective spin susceptibilities, J" jq = 

Aq/(2N) x"f 9 rCt (w ->■ 0), and where N = L/a is the num- 
ber of longitudinal lattice sites (more precisely, x"^'J ct is 
the (a, /?) component of the retarded spin susceptibil- 
ity tensor in the (•) = s, a channel, q denotes the one- 
dimensional momentum along the wire axis and uj is the 
frequency). Like in two-dimensional electron gases^ 
the susceptibilities correspond to the interacting ver- 
sions of the diagrams shown in Fig. [5] (the Luttinger 
liquid formalism however spares us lengthy resummation 
schemes). Similar to the single subband case^ and as we 
will detail below, particle-hole bubbles within each band 
and between the bands result at zero temperature in di- 
vergences of the RKKY exchange couplings at momenta 
2k F i and 2k F 2 in the symmetric channel, and at k F i+k F2 
in the antisymmetric channel (provided that the expo- 
nents defined in the following Eq. ^ fulfill 2gtj < 2), 



Js,i„ - IN X 



cos(2k F iZi) 
ui-21. ri | ±sin(2fc_pi^) 



(10a) 



cos(2k F2 Zi 
+ m 2 k F2 | ±sm(2k F2 Zi 


cos((fc F i + k F2 )zi) 
-fa,*,, =/JVim lF1+t „ ( ±sm((k F1 + k F2 )zi) 





(10b) 



where we used z% = a in. In the semiclassical approx- 
imation employed here, the three magnetizations cor- 
respond to the fraction of microscopic nuclear spins 
participating in either of the three orders. At zero 
temperature, when thermal fluctuations do not weaken 
the order, the magnetizations have to add up to one, 
m 2 k F1 + m 2kF2 + m kF1+kF2 = 1. The helicities and or- 
dering planes are in principle spontaneously and indepen- 
dently chosen for each of the three helices. This analysis 
however neglects the feedback between nuclear spins and 
electrons and is thus not self-consistent. Any finite nu- 
clear spin polarization acts back onto the electrons as a 



static magnetic Overhauser field. Since these fields form 
at the momenta 2kpi, 2kp 2 or kpi + kp 2l they allow 
for backscattering within one band or between the two 
bands. Due to the different inversion symmetries of the 
electronic wave functions in the two bands, however, in- 
traband (interband) backscattering is only possible in the 
symmetric (antisymmetric) channel. The backscattering 
in turn opens up gaps in the electronic spectrum (see 
Sec. lIVB) . If however the electrons are gapped, they can- 
not mediate the RKKY interaction needed to establish 
the nuclear spin order in the first place. 

To find the possible self-consistent orders of nuclear 
spins and electrons, it is useful to note that a helix in the 
symmetric channel at 2kp\ (2fcp2), for concreteness of 
positive helicity, allows only for backscattering between 
right moving spin down and left moving spin up particles 
within the first (second) band. Such a nuclear spin polar- 
ization will thus gap out the right moving spin down and 
left moving spin up modes in the first (second) band, see 
Fig. 0Ja) and Sec. IIVBI The left moving spin down and 
right moving spin up, on the other hand, remain unaf- 
fected. Similarly, we find that an Overhauser field in the 
antisymmetric channel at kpi + &F2 (for concreteness of 
negative helicity) allows for backscattering of right mov- 
ing spin up in both bands and left moving spin down 
in both bands. It can therefore open up gaps in the two 
bands, as depicted in Fig.@|6). These gaps correspond to 
the formation of an electron spin helix in phase with the 
nuclear spin helix (or out of phase for antiferromagnetic 
hyperfine coupling). The Overhauser field induced gaps 
in the electronic spectrum are relevant in the renormal- 
ization group (RG) sense, again provided that the asso- 
ciated exponent fulfills 2gij < 2 (see Sec. lIVB)) . Since g^ 
is intimately related to the Luttinger liquid parameters, 
interactions are crucial for this strong renormalization of 
the Overhauser gaps. 

In the case of a single helix (be it at 2kp\, 2kF2 or 
kpi + kp2), the remaining ungapped electronic modes 
provide the residual RKKY interaction needed to self- 
consistently establish the nuclear spin order. For a su- 
perposition of helices, one finds by counting of the gapless 
right and left moving modes that a self-consistent helical 
order can only occur if either helices at 2kpi and 2kp 2 
coexist (or if there is just one of them) while there is no 
helix in the antisymmetric channel, or if there is a helix 
in this latter channel at momentum kp i + kp2 and no 
helices in the symmetric channel. As an example, con- 
sider the effectively forbidden coexistence of a helix of 
positive helicity in the symmetric channel at 2k p\ and a 
helix of negative helicity in the antisymmetric channel at 
kpi + kpi (see Fig. [4}. This combination of Overhauser 
fields would gap out the entire spectrum in the first band, 
plus parts of the spectrum in the second band. Since the 
first band is fully gapped, scattering events involving this 
band are strongly suppressed. This would eliminate the 
minima at ±2kp i and at ±(kpi + kp 2 ) in the RKKY 
interactions and thus render the analysis inconsistent. 
In addition, the nuclear spins would then reorient into 
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FIG. 4: Opening of gaps around the chemical potential fi in 
the dispersion E(q) due to the magnetic Overhauser field in 
a two-subband quantum wire. Panel (a) depicts the effect 
of a helical field of positive helicity and momentum 2k fi in 
the symmetric channel, which allows for scattering between 
right moving spin down particles and left moving spin up 
particles. The scattering gaps out these two modes. Panel 
(b) shows how a helical field in the antisymmetric channel of 
negative helicity and momentum fe_pi + &F2 gaps out right 
moving spin up particles and left moving spin down particles 
in both bands. As explained in the main text, the intraband 
order depicted in panel (a) and the interband order shown in 
panel (b) are mutually exclusive because their combination 
would gap out the entire lower subband. 



the remaining minima at ±2fcp2, leading to a strong re- 
duction of the magnetic Overhauser fields at 2kpi (sym- 
metric channel) and kp\ + kp 2 (antisymmetric channel) . 
This would finally suppress the electronic gaps due to 
these fields, which is energetically strongly unfavorable. 
The coexistence of a helix of positive helicity in the sym- 
metric channel at 2kpi and a helix of negative helicity in 
the antisymmetric channel at kp\ +kp 2 can therefore be 
excluded. Similar reasonings can be made for other com- 
binations of Overhauser fields. Coexistence of helices of 
equal helicity in the symmetric and antisymmetric chan- 
nel, on the other hand, would lead to competing ordering 
mechanism for the electrons. We thus find that the or- 
ders of the nuclear spins in the symmetric channel (J s ) 
and the antisymmetric channel ( J a ) are mutually exclu- 
sive if the feedback between electrons and nuclear spins 
is taken into account. 



IV. SELF-CONSISTENT SOLUTION WITH 
FEEDBACK EFFECTS 

The ground state of the quantum wire depends on the 
interplay of electrons and nuclear spins, and has to be 
determined as the minimum of the total energy. It should 
thus take into account the energy gain of the nuclear spins 
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FIG. 5: Exponents 2gij of the RKKY exchange couplings, 
see Eq. Q, as a function of vf2/vfi for fixed vfi and in- 
teraction strength U. These exponents describe the regime 
of disordered nuclear spins, e.g. at sufficiently high tempera- 
tures. Consequently, they do not contain any feedback effect 
between electrons and nuclear spins, see Sec. IIVI 



FIG. 6: Exponents 2gij of the RKKY exchange couplings, see 
Eq. ((9}, as a function of vfz/vfi for fixed vfi and interaction 
strength U. On this larger scale logarithmic plot, we also 
show the exponents 2gu,o for decoupled bands, see main text. 
Like for the smaller scale plot in Fig. [5] the nuclear spins are 
assumed to be disordered. 



due to the ordering into the (self-consistently present) 
minima of the RKKY interaction, as well as the energy 
gain of the electrons due to the opening of gaps. 

A. Dominant ordering mechanism: high 
temperature analysis 

To identify the dominant ordering mechanism of the 
nuclear spins (symmetric or antisymmetric channel), we 
first analyze the onset of nuclear spin order upon low- 
ering the temperature. At sufficiently high tempera- 
tures, the nuclear spins are certainly disordered. Using 
\Is,a\/N± ~ /, we infer from the RKKY Hamiltonian 
given in Eq. ([8| that the typical temperature scale for 
the onset of nuclear spin order in the (anti-) symmetric 
channel at momentum q is J s (a),ql 2 - The critical tem- 
perature will be discussed in more detail in Sec. [V] Just 
before the onset of the nuclear spin order, when there is 
no Overhauser field present yet, the RKKY interaction 
mediated by the fully ungapped electron sector can be de- 
rived from the electron Hamiltonian given in Eq. (llb[) . In 
order to include the density-density interactions, we treat 
this Hamiltonian with standard bosonization techniques, 
which allow for its full diagonalization, see Appendix |B| 
This diagonal form finally allows the determination of 
the spin susceptibilities and thus the RKKY exchange 
couplings, as detailed in Appendix [Cl 

The resulting exponents 2 t/jj defined in Eq. (JSJ , which 
control the strength of the RKKY interaction, are plot- 
ted in Figs. [5] and [5] as a function of vpijvFi for fixed 
vfi and fixed interaction strength U (thus following the 
filling of the second subband). The interaction strength 
is chosen such that the Luttingcr liquid parameter in the 
charge sector of the first subband is K& = 0.5, while 
we use K a i = 1 because the interactions in the spin sec- 
tors are assumed to be much smaller than the ones in 



the charge sectors (in agreement with experiments) . 19 i 20 
We find that at the special point vpi = i'F2, all sus- 
ceptibilities have the same exponent. This has already 
been discussed in the context of Carbon nanotubes, 10 
where the two subbands correspond to the two inequiva- 
lent Dirac cones. Away from this special point, however, 
the exchange coupling J s associated with the symmetric 
superposition of nuclear spins I s is always more singular 
than the exchange coupling J a associated with J a (we re- 
call that J s ~ X11+X22 while J a ~ Xi2 + X2i)- Therefore, 
the system will initially order in the symmetric channel. 

Since the spin susceptibilities and thus the RKKY 
exchange couplings are mediated by particle-hole pairs 
exchanged between two nuclear spins, see Fig. [2] it is 
physically not surprising that the interband RKKY ex- 
change can only be effective if the particle and the hole 
propagate at comparable velocities. Consequently, the 
two subbands are essentially decoupled when w^i > vfi 
or tifi <C vp2- This is shown in the larger scale plot 
in Fig. [5J which indicates that the exponents 2gn ap- 
proach their decoupled values 511.0 and 522.0 (calculated 
in the absence of interband density interactions) away 
from vfi ~ VF2- Only close to vfi — VF2, the interband 
channel participates in the RKKY interaction. At the 
same time, the weight of the intraband exchange chan- 
nels decreases and the corresponding exponents gn and 
522 are larger than in the decoupled case. A related sit- 
uation arises in Coulomb drag setup, where the drag is 
most efficient if the Fermi velocities of the two wires are 
approximately equals 

B. Self-consistent low temperature ground state 

Given that the dominant instability at high temper- 
atures occurs in the symmetric sector, we calculate the 
low temperature ground state assuming that the nuclear 
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spins form helices in the symmetric channel. This as- 
sumption will be justified a posteriori in Sec. [VT] We 
start from a general superposition of nuclear spin helices 
in the symmetric channel, 



(11) 







f cos(2k F iZi)\ 


INj_ 


m 2 k F1 
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with magnetizations m,2k F1 + m2k F2 = 1 at zero temper- 
ature, and subsequently minimize the total energy with 
respect to the magnetizations rri2k F1 and m2k F2 - The 
system is analyzed in a bosonized language, in which the 
annihilation operator for a particle of spin a in band i 
reads 



r{z) 



-W) 



V2t. 



(12) 



where r — R, L = +, — , the corresponding Klein fac- 
tors are denoted as t/ r j CT , and a is a short distance cutoff 
(here taken to be the lattice spacing) . The bosonic fields 
4>i a and 9i a fulfill the standard commutation relation 
[(/)ia(z),8i'a'(z')] = ki'&oo> (w/ 2 ) sgn(z' - z). 27 It is fur- 
thermore convenient to introduce spin and charge degrees 
of freedom via the canonical transformation <j>. c (z) — 



9, t 



± <Pn)/V2 and d.c (z) = (6 n ± 8 tl )/V2. 



As shown in Appendix [D] the hyperfine coupling be- 
tween the ordered nuclear spins and the electrons yields 
various cosine-terms. If the Fermi momenta fc^i and kp2 
are non-commensurate, the only non-oscillating terms 
read 



H. 



non— osc. 



E 



B 



xy,i 



2na 



dz cos 

ic \ "i 



■>) 



(13) 



where B xy ^ = I Aom2k Fi - These Overhauser fields are 
RG relevant and thus lead to an ordering of 4>i C + 0j s , 
i.e. half of the electronic degrees of freedom. This in 
turn feeds back to the nuclear spins by a modification of 
the RKKY exchange interaction, which is now only me- 
diated by the gapless part of the electronic spectrum. To 
include this feedback effect into the theory, we start from 
the non-diagonal version of the electron Hamiltonian in 
the presence of the Overhauser fields and interband in- 
teractions, 



*e=E E 

2—1,2 j — c,s 

dz AU 
2n 7r 

Bxy,i 

2ira 



dz ( u 



2tt \K 



{d z (t>ij? + u l3 Kij {d z 6 v] )'- 



E 



(dAic)(dM (14) 
dz cos (y2((j) ic + 9 is y\ . 



Here, Uij and Kij are the effective velocities and Lut- 
tinger liquid parameters in the spin and charge sectors of 
the two bands with Uij = vpi/Kij. This representation 
indicates that the interband Coulomb interaction and the 
intraband Overhauser fields have competing effects and 
should be treated simultaneously. While the interband 
interaction is diagonalized by band-mixed boson fields 
given in Eq. (|B4[) . the magnetic Overhauser fields induce 
an intraband ordering and gap out parts of the spectrum 
within each band, which opposes a complete mixing of 
the bosonic fields. Fig. [7] depicts this in terms of RG 
fixed points. The first line of Eq. (|14p corresponds to 
the fixed point a with two decoupled Luttinger liquids. 
Adding an interband interaction U\2 — U but no cou- 
pling to the nuclear spins takes the system to a different 
fixed point b described by 2 decoupled Luttinger liquids 
with modified Luttinger parameters and velocities whose 
Hamiltonian is given in Eq. (1B7[) . This flow is shown in 
the upper panel of Fig. [7] The fixed point b has been 
used to calculate the spin susceptibilities in Sec. IIV Al 
The lower panel of Fig. [Jj on the other hand, depicts how 
the RG flow is changed due to the hyperfine coupling Aq 
to the nuclear spins. Through the Overhauser fields, the 
hyperfine coupling gaps out half of the electronic degrees 
of freedom. As a consequence, the fixed points a and 
b are not stable at low temperatures. If there were no 
interband interaction, the system would be described by 
a fixed point c with 2 gapless and 2 gapped Luttinger 
liquids. Now adding an interband coupling finally takes 
the system to yet another fixed point d, which is the one 
we are after. 

Given that the intraband Overhauser fields are 
strongly RG relevant, they dominate the initial RG flow 
compared to the marginal interband density-density in- 
teraction. The latter is in any case only important when 
Vpx ~ vf2 and may then indeed slightly modify the scal- 
ing of the Overhauser gap, but not the presence of an 
ordered state. We thus first describe the flow from fixed 
point a to fixed point c and subsequently take into ac- 
count the interband interaction. Close the to fixed point 
a, the term cos(v / 2(</>ic + Sis)) has the scaling dimension 
cji = (Ki C + X/Ki s )/2, Parameterizing the running short 
distance cutoff as a(b) 
the fixed points a and c reads 
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d\og(b) 



0--9i)B. 
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FIG. 7: RG flow of the two-subband quantum wire system. 
The upper panel (a) shows the flow for density-density cou- 
pled subbands in the absence of hyperfine coupling. For any 
value of the interband interaction U12 , the system can be de- 
scribed by 2 decoupled spinful Luttinger liquids (LL), charac- 
terized however by different Luttinger liquid parameters and 
velocities. The lower panel (b) shows how the flow is modified 
due to the presence of the hyperfine coupling and a result- 
ing ordering of the nuclear spins in the symmetric channel. 
Here, the interband interaction would correspond to the out- 
of-plane axis (not shown), and panel (6) is the projection of 
the full flow into the given plane. 



criterion, reads 



KiB 



xy.i 



u i+ 



KiaB x 



(l-fli)/(2- ff< ) 



(18) 



If however the now stops at b* L T = a* /a with a{b* L T ) = 
a* = min{L, Ui + /T} due to the finite system size or 
temperature, the expansion of the cosine to second or- 
der yields a gap smaller than the finite size gap (or the 
temperature) m+/a*, and fluctuations beyond Gaussian 
order may be important for a given physical observ- 
able. Nevertheless using the Gaussian approximation, 
the Overhauser field induced gap can be expressed as 




Kia(b* L T ) B xy ,i(b* L T ) 



(19) 



The renormalized Overhauser fields at the end of the 
flow, B* yi = B XVi i(viim{b* ,b* L T }), can be interpreted 
as renormalized hyperfine coupling A* = Aq B* y i /B xy ^ 



for the 



IJ i + 



modes (we recall that B xy 



A ). This 



renormalized hyperfine coupling expresses how the or- 
dered electronic modes collectively couple to the the po- 
larized nuclear spins. 



C. Residual nuclear spin Hamiltonian 



As discussed in Appendix [E] (and similar to the single 
subband case)^ this leads to gaps for the bosonic 
fields </>i+ ~ fa c + 9i s . There are however two more 
bosonic modes fa- resulting from the canonical trans- 
formation 4>c,9s —> <P+,fa-- These modes do not couple 
to the Overhauser fields and remain gapless. 

Expanding the cosine-terms to second order— and us- 
ing u i+ = (u ic K ic + u is /K is )/(K tc + 1/JQ S ) and Ki = 
2gi = K ic + 1/K is yields 



A i+ (6) = 



' KiU i+ B xy: i(b) 



a(b) 



(16) 



Assuming the presence of two nuclear spin helices of 
the form (Jl]) at low temperatures, we found that the 
fields 4>i+ ~ 0ic + Ois and (j>2+ ~ fac + #2s are gapped. 
Rather than through an exchange interaction mediated 
by all electronic modes, the nuclear spins feel the pres- 
ence of the electrons through an RKKY interaction me- 
diated only by the remaining gapless fa- modes, as well 
as through the finite magnetization of the electron spins 
of the 4>i+ modes, which form helical spin density waves 
locked to the nuclear spin helices^ To make our cal- 
culation self-consistent, we thus have to recalculate the 
RKKY interaction mediated only by the fa- modes. As 
detailed in Appendix [F] we obtain the residual exchange 
couplings in x and y direction as 



The flow in band i stops if either the running gap reaches 
the running energy cutoff, A i+ (b) = Ui+/a(b), or if the 
short distance cutoff equals the wire length, a(b) = L, 
or if the finite temperature cuts off the divergences at 
a(b) = u + i/T. The first criterion defines the RG stage 
b* as 1 = Ki a(b*) B xy j(b*)/m + , and thus 



jxxl _ jyyl _ jl , 
s,q s,q s,g,l _r 

jxxl _ jyyl _ jl 

u a.,q — u a,q ~ ° a,q,12 



J s,q,2 ' 
1" J&,q,2l ■ 



(20a) 
(20b) 



At finite temperatures, the RKKY exchange couplings in 
the symmetric sector read 



b* = 



KiaB x 



1/(2-3.) 



(17) 



The physical gap A* + = Ai + (b*) in the electronic spec- 
trum, provided that the RG flow is stopped by this first 
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FIG. 8: Exponents 2g' i j of the residual RKKY exchange cou- 
plings after the formation of nuclear spin helices at momenta 
2£;.fi and 2kFi, as a function of vf2/vfi for fixed vfi and 
interaction strength U such that K\ c = 0.5. 



where T is the standard Gamma function, Uj_ = 
(uic/Kis + Uis Ki e )/(Ki c + 1/Kis) and where the expo- 
nents are given in Fig. [5J The intraband couplings 
have analogous expressions. At zero temperature, this 
translates to the divergent behavior 



■CU» F1 ~l«=F2JI*i 



jxxf I 

J s-q \ q ai±2k F2 



|g=F2£;F- |2 ^ 2 ~ 2 



TM/ I 

"^a.g l g «±(fe F1 +fc F2 ) 



«F1| 

|<7T (&fi 



fc F2 )l 29 -- 2 
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Note that for the calculation of the new exponents g[p 
we have neglected renormalizations of the Luttinger liq- 
uid parameters and velocities within each band due to 
interband interactions during the flow between the fixed 
points a and c in Fig. [7] Since this flow is rather short 
(we recall that the Overhauser fields are strongly RG 
relevant), and since interband interactions have a small 
effect unless vfi ~ vfi, we expect this approximation 
not to change our results. 

Like for the single band wire, the RKKY interaction 
is now anisotropic. Correcting an error in Ref. [Toj . the 
RKKY coupling along z remains non-singular in the he- 
lically ordered state. This is an can be understood from 
Fig. 21 To obtain a divergent contribution to the z com- 
ponent would require backscattering from a left-mover to 
a right-mover without spin flip. Such a process is how- 
ever not possible in the helically ordered state where the 
right moving spin down and the left moving spin up par- 
ticles are gapped. Technically, the spin susceptibility in 
z direction is regularized by the fact that it involves the 
field 6*+ canonically conjugate to the ordered field <f>+, 
which cuts off its divergence ! 13 ' 28 " — 

To conclude, we find that the effective Hamiltonian 
for the nuclear spins including the feedback between the 
latter and the electrons reads 



^RKKY 



jololI J eta. I 

TCt S ?ff TLX I TLX S i Q TOl. 

' G ~ q N 2 ,q a '~3 ]\J-2 J a,q 



(23) 



(A\ (S lhlll ) + A* 2 (S 22 , H ))-I S , H 



N 



As has been discussed in the single subband case, the 



effective magnetic fields 



acting on the nuclear 



spins due to the polarization of the (partially) ordered 
electron system are however negligible compared to the 
residual RKKY interaction and will therefore be omitted 
in the following^ 



V. COEXISTENCE OF HELICES AT 2k F i AND 
2k F2 AND ONSET OF NUCLEAR SPIN ORDER 
AT FINITE TEMPERATURE 

Having established the effective theory in the pres- 
ence of both hyperfine coupling and interband interac- 
tions, we are finally in the position to analyze the ground 
state of the system at low temperatures, where the nu- 
clear spins order in a superposition of two helices at mo- 
menta 2kpi and 2kp 2 with magnetizations m 2 k F1 and 
ni 2 k F2 = 1 — rn 2 k F1 at zero temperature, see Eq. (jlll) . 
The coexistence of these two helices can be inferred from 
the fact that the energy loss of orienting parts of the 
nuclear spins into the shallower minimum of the RKKY 
interaction, see Fig. [3l is overcompensated by the en- 
ergy gain due to gapping out a second electronic sector. 
Considering for example kp 2 < kpi, the deepest mini- 
mum of the RKKY exchange is at momentum 2kp 2 . The 
loss in both nuclear spin energy and the electronic gap 
of the second band while populating the second mini- 
mum with a small magnetization m 2 j CF1 scales linearly in 
Ln 2 k F1 , while the gain in energy due to the gap in the first 

subband scales as w^ 2 > m 2 k F1 - A similar argu- 
ment does not hold for the helical interband polarization 
Lnk F1 +k F2 because the latter would gap out the electronic 
spectrum and suppress the RKKY interaction altogether, 
see Sec. IIII1 The ground state magnetizations can then 
be obtained by minimizing the total energy with respect 
to m 2 k F11 which depends on the ratio of vp 2 to vfi, as 
well as the temperature (through the RKKY interaction). 
The nuclear spin order is also stable if interband spin 
density interactions or current-current density interac- 
tions were non-negligible. Being quadratic in the bosonic 
fields, these interactions could be taken into account by 
a generalized transformation diagonalizing the interband 
interaction after the helix formation in Sec. IIV CI and 
would therefore not modify the presence or absence of 
electronic gaps. 

At larger temperatures, thermal fluctuations decrease 
the nuclear spin polarizations m 2 k F1 and m 2 k F2 , which 
also has to be fed back into the minimization procedure. 



10 



When the temperature is eventually high enough, ther- 
mal fluctuations entirely suppress the nuclear spin order. 
Very roughly, this is expected to occur when the thermal 
energy is larger than the RKKY exchange at the helix 
momentum. An upper bound to the critical tempera- 
ture is defined by the residual RKKY interactions that 
profit from the reduced exponents g'^ defined in the or- 
dered phase. This yields the estimate fcsT C j ~ I 2 J*f 
for the transition temperature T C i in band i. Since the 
two helices at momenta 2kp i and 2kp2 are sustained by 
different minima of the RKKY exchange coupling with 
different depths, two distinct ordering temperatures can 
be observed. We derive these ordering temperatures in an 
approximation known as the Tyablikov decouplingi 31 i 32 
This method allows the self-consistent determination of 
the critical temperatures associated with the two helical 
nuclear spin orders in a mean field type approach. We 
find that the helical order in band i sets in for tempera- 
tures below 



dq — 



a-i{q) 



with 



o)- Mq)- JIM] 



(24) 



(25a) 
(25b) 



with Ai(q) = J*f (g + 2k Fi )/2 + J™'{q - 2k Fl )/2 and 
Jsi(q) being the contributions of band i to the total spin 
susceptibility (see Eq. (|20ajl ). Details of the derivation of 
the critical temperatures T C i can be found in AppendixlGl 
The spin susceptibilities are dipped around q — ±2kpi, 
see Fig. [3j and we thus have a 2 — f3 2 ~ q 2 close to q = 
0. The integral in Eq. (|24|) is therefore infrared (IR) 
divergent. This divergence is cut off by the finite size of 
the sample at the momentum \q\ ~ 2ir/L. For not too 
long quantum wires, the integral can be approximated 
by the implicit equation (note that J^f has a power law 
temperature dependence) 



T 



(26) 



which confirms our order of magnitude guess. These 
critical temperatures T C i found with the Tyablikov de- 
coupling are consistent with the result obtained in a 
magnon analysisi 10 i 33 For longer wires, however, one 
needs to compare the critical temperature given in 
Eq. (|26p . resulting from the non-singular part of the 
integral in Eq. (|2"4")l . to the IR divergent contribution 
~ [(d 2 a/dq 2 )di P ]^ 1 La/ir 2 . For L -» oo, the diver- 
gent contribution dominates, and the critical tempera- 
ture goes to zero. This limiting behavior is consistent 
with the general theorem that there is no order possible 
in these systems at finite temperatures in the thermody- 
namic limit i^i 



We thus find that the nuclear spins in a finite size quan- 
tum wire with two subbands form two distinct helices 
of momenta 2kpi and 2kp2, and that these two orders 
have different critical temperatures. The superposition 
of helices with different pitch lengths results in a beat- 
ing pattern for the nuclear spin polarization. Since each 
helix order gaps out one electronic mode, the quantum 
wire exhibits a two step reduction of the conductance 
from 4e 2 /h to (approximately) 3e 2 /h at the first order- 
ing temperature, and to (approximately) 2e 2 /h at the 
second ordering, which provides an experimental signa- 
ture of the double helix order. 



VI. INTRA BAND VS. INTERBAND ORDER 

A number of consistency checks of the results obtained 
here can be directly translated from the single band case, 
and we refer the reader to Ref. .10] for details. What re- 
quires a little more care is the stability of the order in 
the symmetric channel with respect the order in the an- 
tisymmetric channel. The assumption of an ordering in 
the symmetric channel was based on the more divergent 
spin susceptibility of the latter with respect to the disor- 
dered fixed point b in Fig. [JJ or alternatively the higher 
RG relevance of the Overhauser field associated with this 
channel close to this fixed point (note that close to the 
fixed point a of Fig. interband scattering processes are 
even irrelevant in the RG sense). If the ratio of Fermi 
velocities becomes of order one, however, the two orders 
have comparable scaling dimensions. One can suspect 
that in this regime, also the electronic gaps due to either 
possible ordering, as well as the nuclear spin energies (set 
by the RKKY interactions and thus also controlled by the 
exponents gij that determine the scaling dimensions of 
the Overhauser fields) will become of the same order, and 
that the ground state can only be inferred from a detailed 
comparison of the various possible nuclear spin orders 
and associated electronic states. Similar considerations 
apply to the limit of weakly interacting electrons, where 
all gtj — > 1. For smaller vp2jvp\ and experimentally ob- 
served strongly interacting Luttinger liquids, however, a 
helical order of the nuclear spins in the antisymmetric 
channel can be excluded. 

For vp2 < 0.5 Vpi, specifically, the stronger residual 
spin susceptibility in the symmetric channel, see Fig. [51 
ensures a helical ordering at 2k p\ and 2kp2 at low enough 
temperatures. For larger vp2/vp\, our analysis becomes 
inconsistent for T — > when the divergent spin suscep- 
tibilities overrule any finite gap in the electronic sectors. 
Given however that the maximal ratio of Fermi veloci- 
ties before a third band is filled is Vp^/vpi ~ 0.7 (for 
two subbands in a single wire) , most of the two subband 
regime will exhibit the formation of two helices at 2kp\ 
and 2kp2, while the regime vp2 — > vp\ deserves further 
analysis elsewhere. 
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VII. SUMMARY 

In this work, we have considered the interplay of elec- 
trons and nuclear spins in a two-subband quantum wire. 
Similar to a single subband wire, the hyperfine coupling 
between electrons and nuclear spins leads to an ordering 
of the nuclear spins and consequently a gap for half of 
the electronic spectrum~*iP This gap is strongly renor- 
malized by interaction effects. Different from the single 
subband case, the nuclear spins in a two-subband quan- 
tum wire form a superposition of two helices with distinct 
pitches Ai = ir/kFi and A2 = Tr/kF2, which gives rise to 
a beating pattern in real space. The two helical orders 
set in at different temperatures. As a result, the con- 
ductance of the wire exhibits a stepwise reduction from 
4 e 2 /h at highest temperatures to (approximately) 3 e 2 /h 
after the formation of the first helix, and finally to (ap- 
proximately) 2e 2 /h at the formation of the second helix. 
This provides an experimental signature of the double 
helical order in two-subband quantum wires. 

To establish these results, we have analyzed the quan- 
tum wire in a self-consistent analysis based on the well- 
separated dynamics of electrons and nuclear spins. As 
has been shown in SecHH the hyperfine coupling between 
electrons and nuclear spins can be decomposed into a 
symmetric and an antisymmetric channel. The electrons 
thus mediate two kinds of RKKY interactions between 
the nuclear spins. As has been discussed in Sec. lIIIl these 
different interaction channels compete for the formation 
of nuclear spin order. The final ground state has then 
been derived in a self-consistent approach that takes into 
account the interplay of nuclear spins and electrons by 
describing the nuclear spins as an effectively static Over- 
hauser field for the electrons, while the latter mediate an 
effectively instantaneous RKKY interaction for the for- 
mer. This calculation has been carried out in Sec. IIVI 
We find that the two-subband quantum wire preferably 
orders in the symmetric hyperfine channel, which leads 
to the formation of the two superimposed nuclear spin 
helices of pitches Ai and A2 at lowest temperatures. At 
higher temperatures, each helix is associated with its own 
critical temperature, implying the above described reduc- 
tion of the conductance in two steps as the temperature 
is lowered, see Sec. IVl 
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Appendix A: RKKY interaction 

In order to derive the RKKY interaction by integrating 
out the fermionic degrees of freedom, we first rewrite the 
hyperfine coupling in momentum space and imaginary 



time. Using (r) = Y^. e~ iqZi I^ )q (t) with 



(0« v 



where in denotes the longitudinal nuclear spin lattice site 
index and with a being the lattice constant, we obtain the 
action corresponding to H ne as 

Snc = TF-Jr / dT^2(S n ,q( T ) + S 22,q( T )) ' I s,q( T ) 

(AI) 

+ ^ / drJ2 (Si2, q (r) + S 21 , q (r)) • I a , ? (r) , 
± q 

where N = I/a is the number of lattice sites 
in longitudinal direction, and with Sij, q (r) — 

Ea,/9,fc c Lfc+g( T )( <T ^/ 2 ) c i«fe(' r )_- The RKKY interac- 
tion can be obtained by expanding the total action S — 
S c + 6> ne to second order in S ne and integrating over the 
fermions. Re-exponentiation then yields the action 



5rkky 



4lJ_ 

Nl 2N 2 



(It I tlr'Y, 



(A2) 



x (i:-q(r)xt P q(r-r')I^ q (r') 
+ ^-q(r)x^q(r-r')I^ q (r')) 
with the imaginary time spin susceptibilities 

Xf q (r - t') = (S? 1 ,_,(r)Sf 1 ,,(/) + S^_ q (r)S^ q (r')) , 

(A3) 

<(t - r<) = {S^_ q {r)S^ q {r') + S^_ q {r)S( %q {r')) . 

(A4) 

The averages are taken with respect to S e and evaluated 
in Appendix [Cj Given that the nuclear spin dynamics are 
associated with much longer time scales than the electron 
dynamics, we can approximate (with j3 = 1/T in units of 
k B = 1) 

X(-u(r -r') = - Y J e -^ (T - r,) X ( .), g K) (A5) 

= S(t - r')X(.),g(w n = 0) . 

The real time Hamiltonian corresponding to the action 
5rkky can finally be found by analytical continuation to 
the retarded spin susceptibilities, 



A 2 1 

- 1 - q 

X (l?,-q(t)xf q let ("^0)ie q (t) 

+ 4 a - 9 Wx^' rot (^^o)/f i9 w) 



(A6) 
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Appendix B: Diagonalization of the electron sector 



In order to diagonalize the electronic Hamiltonian H e , 
we first bosonize the annihilation operator for a particle 
of spin a as 



where 



r{z) 



R,L 



U, 



V2t 



•'<"' e -i{r<t>j<r(z)-ej„(z)) 



(Bl) 



-, — , the corresponding Klein fac- 



tors are denoted as U r j a , and a is a short distance cut- 



off (here taken to be the lattice spacing). The bosonic 
fields 4>j a and 6j a fulfill the standard commutation re- 
lation [<j>jo{z),0ji a >(z!)] = 5jj> 5aa' (nr/2)sgmV — z) , 
Since the interaction U in Eq. ([T]) is of density-density 
type, it can be taken into account by an exact basis trans- 
formation. Because furthermore the Fermi velocities vpi 
and Vf2 in the two bands are different, it is most con- 
venient to first diagonalize each subband separately by 
introducing spin and charge degrees of freedom via the 
canonical transformation <p.c(z) — {4>*\ ± 0,^)/v / 2 and 

^^(z) — (9ft ± 6ft)/*/2. In this basis, the electronic 
Hamiltonian reads 



^ = E E 

i=l,2 j— c,s 

dz AU 

2lT 7T 



£ fe(^-) 2 +««^i(w 2 



(d z 4> lc )(d z (t>2c), 



(B2) 



where the effective velocities and Luttinger parameters 
are as usual given by 




Kir. = . . 

irvFi 

K is = 1 , u is = v Fl 



v Fi fJ30 \ 

— , (B3a) 
(B3b) 



The diagonal electronic Hamiltonian is now obtained by 
a second canonical transformation 



He = 



VFl 



u c+ (l + Al) 




A\ VFl 



I A 2 V F 2 

r u c+ 

V F l(l + A 2 c ) 

I A 2 C u c+ 

, VF2(1+A 2 ) 

with the velocities 



72c 




u c -(l + Al) 



VF2 

u c .(l + Al) 



A 2 C u c - 
vfi(1 + A 2 c ) 



2\ Vc ~ ' 



(B4a) 



(B4b) 
, (B4c) 



7 e+ 



VF2(l+Al) 



Be- ■ (B4d) 



U c ± = 



\ 



u 2 + uL 



± 



2U 



VF1VF2 



(B5) 



and with 
A c 



(4U/n) ^Jv F ivf2 



\J ( u ic - u lcf + (4U/ir) 2 VpiVF2 + uj c - u\ c 

(B6) 



27 The electronic Hamiltonian can then be written as 



fc=± 

E 

.3=1,2 



H ° = E If / dz {( d ^kf + (d z e ck f) (137) 



VFj 

2tt 



dz {^d z 4> sj f + (d z 6 S3 y 



We note that a generalized form of the transformation 
given in Eq. (IB4[) would also allow to take into account 
the spin density-density interaction, charge current- 
current interaction and spin current-current interaction 
neglected here. 



Appendix C: Evaluation of the spin susceptibilities 

The real time RKKY Hamiltonian given in Eq. ([8]) de- 
pends on the zero frequency components of the retarded 
spin susceptibilities. The latter can most conveniently be 
derived starting from the imaginary time susceptibilities 
given in Eqs. (|A3j) and (|A4|) . We obtain 



<sW^>')> = E E f^r 

k,k r a,b,a',b' 

X (4ak- q ( T ) C jbk( T ) C l'a'k'+ q ( T ') C fb'k'( T ')) 
J 

E 

,b,a' ,b 

j 



(Cl) 



Zj±Z^L L j d(z-z')e lq{ " ' 



x{cUz,T)c jb {z,T)clA^T% lb ,(z'y)) . 

Since the Hamiltonian conserves the total spin, these av- 
erages (evaluated with respect to H e ) are diagonal in spin 
space. In addition, we find that forward scattering con- 
tributions have a non-singular spin susceptibility, while 
backscattering events result in divergences. We therefore 
only evaluate the backscattering terms. To this end, the 
electronic Hamiltonian is bosonized and (depending on 
the fixed point one is interested in) diagonalized as de- 
tailed in Appendix [EJ The divergent part of the x and 
y components of the spin susceptibility within the first 
band are identical and read 
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(C2) 



4(27ra) 2 



d(z - z')e iq{z - z ' ] { e - i2kFl( - z - z ' ) 



1^2. The z component finally has the expression 

x \ 2 q = E T / d ( z - z')e tt - q - r(kF1+kF2 ^ z - z ^ (C5) 

x (r\ a (z, r)f 2(T (z, r)f| CT (z', r> 1(T (z', r')) 

1 d(z - z V 9 ^"^ ^ ,(tF1+fcF2)(z " z,) 
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x e 
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with r = R, L — + 
z-component reads 



The backscattering part of the 



Xn, 



= Ei /*(*-*')« 



(C3) 



x (rj^z, r)r 1(T (z, r)r f 1(7 (z', r> 1(T (z', r')) 
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L 
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+ ( e i ^^ lc ( 2 ' T ' _ ^ ls ( Z:T )~^ lc ^'' r ') + ^ ls ( z '' T '' ) ' ) ) +hc j 

The spin susceptibility within the second band can be 
obtained from the one in the first band by exchanging 
the index 1 — > 2. The spin susceptibility between the 
bands finally has the divergent part 
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In addition, there is a contribution xfi 9 = X21 g which 
can be obtained from X12 g = X12 9 by the index swap 



x e 



x e 



i(l/v / 2)(e lc (z' ,r')+0i a (z' ,T')-e 2c (z' ,T')-e 2 »(z' ,t'))\ 
a i(l/v / 2)(0lc(z,T)-0 ls (z,r)+0 2c (z,T)-02 3 (z,r)) 
-i(l/v / 2)(ei c (z,r)-e l3 (z,r)-e 2c (z,r)+e2 s (z,T)) 
-j(l/v / 2)(0lc(z',r')-0i s (z',r')+0 2 c(z',r')-<!i2 s (z',r')) 



x e i(l/v / 2)(ei c (z',T')-ei s (z',r')-e 2c (z',r')+92 ii (z',T'))^ 

+ (kFi -> -k F i,4>ij -> -4>tj)j , 

plus again an additional contribution that is obtained 
under 1 O 2. These averages are now evaluated with 
respect to the diagonalized Hamiltonian H e given in 
Eq. (|B7[) . We explicitly state the calculation of xfn the 
other susceptibilities can be evaluated analogously. Fol- 
lowing Ref. and using Eq. (|B4al) . we find that 



X ^ q = -j2^r ' d{z ~ z ' ] 



(C6) 



^ e »(5+2fepi)(*-«') _|_ e i(q-2k F1 )(z-z')^ 



x e 
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A = 



VFl 



u c+ (l + A 2 c ) 



B=o A l VF \^ , (07) 



ti c -(l + .A; 



21 ' 



and where the function i* 1 ^ (z, r) is at zero temperature 
given by 



F 1 ( ' ) (z,r) = log 



f > /r' + (u(.)lr| + a) ii 



(08) 



where is the velocity associated with the respective 
mode. At finite temperatures, it becomes 



F«(^r)=log(^«,/sinh 2/ 



a I rSln \ a 
PU(.)J VP 



(09) 
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(where /3 = 1/T in units of ks = 1). Being mainly inter- 
ested in the scaling behavior of susceptibilities, i.e. their 
power law exponents, one can neglect the difference in 
velocities and obtains 



Xn, 



(27ra) 2 



ira/ (0u) 



sinh' 



sin 



2 S tt(t-t') 



(CIO) 



Here, u denotes the common velocity of the different 
modes in this approximation. The physical response of 
the system can now be obtained by analytic continuation 
of Eq. ([CIO]) . The latter is described in detail in Ref. [27| 
and yields 



a:, ret I 



2-2 ff i 

47r 2 u \ 2ira J 

r(^-i^(9 + «2* F1 )) 
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Er 2 (i- Su ) 



(Cll) 

2 



where T is the standard Gamma function and with 



5ii 



A 2 + B 2 + l/K u 



(C12) 



At zero temperature, this translates into 



a;, ret 

Xii 



1 lw-»0 



L sin(7rgn) 

47T 2 W 



E 

K=± 



2-2 S n 



(C13) 



The spin susceptibility within the first band xn thus di- 
verges with a power law exponent 2— 2gn at the backscat- 
tering wave vector 2kp\. We find similar power law diver- 
gences for X22, X21, and X12 characterized by exponents 
2 — 2gij. These exponents are plotted in Figs. [5] and |6] 
Due to the transformation detailed in Eq. (|B4I) , they de- 
pend on the interaction strength between the two bands 
as well as the ratio of Fermi velocities vp\/vp 2 . 



Appendix D: Overhauser fields and their coupling to 
the electronic modes 



In this Appendix, we shortly review how an Over- 
hauser field of one of the two effective nuclear spins af- 
fects the electrons in the quantum wire. Let us first start 
with a helical field of positive helicity in the symmetric 
channel at 2k p\, 



/ 'cos(2kpiz) 

Bov, s (z) = (I s . z ) = B I sm(2kpiz) 





(Dl) 



This Overhauser field is coupled to the electrons by the 
hyperfine interaction, 



H en = p- Jdz (S n (z) + S 22 (z)) ■ B 0v , s (z) (D2) 



AqB q 

2N 



e - i2kFlz c u + h.c.) 



'2-+. 6 " ^24, I li.C. 



The only non-oscillating cosine-terms deriving from this 
interaction (unless the two Fermi wave vectors are com- 
mensurate) are backscattering terms of a right-moving 
spin down and a left moving spin up within the first band. 
Dropping the Klein factors that cancel out when taking 
averages subsequently, the bosonized version of this term 
reads 



H, 



non— osc. 



AqBq 

2-KaN± 



dz cos 



(V2(4>u 



+ 1 



(D3) 



This term is strongly RG relevant and gaps out the com- 
bination <p c i + 8 s i in the first band, corresponding to 
right moving spin down particles and left moving spin 
up particles (see Ref. [To| for a detailed discussion). An 
Overhauser field of negative helicity would gap out the 
combination <fi c i — 9 S \ that corresponds to right moving 
spin down particles and left moving spin up particles. 
A helical Overhauser field with momentum 2k p 2 , on the 
other hand, would gap out the corresponding combina- 
tions <j) C 2 ± 9 C 2 in the second band. 

We now turn to an Overhauser field in the antisym- 
metric channel. For 



cos((k F i + k F2 )z) 
Bov,a(z) = (I^z) = B ( -sin((fc F i + k f2 )z) ] , 

' ,D4) 



the hyperfine interaction is 



H 



Ao_ 
A Bo 



J dz (S 12 (z) + S 21 (z))-B 0v ^z) (D5) 
2N± J dz (^e^+^c^ + h.c.) 

+ (4t^ 1+fcF2) ^u + h-c-) ■ 

Bosonizing this Hamiltonian, we find two non-oscillating 
cosine-terms. The first line allows for backscattering be- 
tween right moving spin up particles in the second band 
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and left moving spin down particles in the first band, 
while the second one allows for scattering between right 
moving spin up particles in the first band and left mov- 
ing spin down particles in the second band. Bosonization 
finally yields 



rrnon-osc. AqBq 

2iraN 



AqBq 

2iraN± 



where 



dz COS fV2(012c - 0l2s)) 
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(E2a) 
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(E2c) 

(E2d) 



Kuy/Ki'" V K is Ki 
with Ki = Ki c + 1/Ki s = 2gi. This transformation yields 



<\>A + <t>j\. 
V2 



V2 



(D7) 



Like the cosine-terms within a single band, the interband 
cosine-terms are RG-relevant if the exponent 2gi% defined 
in Eq. is smaller than 2 and open up gaps for 4>i2 C — 
9\2 S and </>2ic — #2is, or equivalently for right moving spin 
up particles in the both bands and left moving spin down 
particles in both bands. A field of opposite helicity would 
have gapped out 4>i2c+8i2s and 4>2ic+&2is, corresponding 
to left moving spin up particles in both bands and right 
moving spin down particles in both bands. 



Appendix E: RG flow due to the Overhauser fields 

Taking into account the Overhauser fields, the elec- 
trons arc subject to the effective Hamiltonian 



2=1,2 j—c..s % 3 

J 2tt 7T 



(El) 
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i=l,2 
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27ra 



dz cos 



As discussed in Sec. IIVB1 the Overhauser fields B xy>i 
and the interband interaction U have competing effects 
and, when taken on their own, drive the system to two 
different fixed points. Their combined effect is taken into 
account by first analyzing the flow due to the strongly rel- 
evant Overhauser field whilst neglecting the flow due to 
the marginal interband interactions, which will be treated 
in a second step. 

At the initial fixed point described by the first line of 
Eq. (|E1[) . the cosine-terms of the Overhauser fields have 
the scaling dimensions ^ = (Ki C +l/Ki S )/2. To integrate 
the flow of B X y t i, we perform a canonical transformation 
that switches from spin and change degrees of freedom 
to fields 4>i+ ~ 4>i c + Bis whilst preserving the scaling 
dimension of the cosine-term with respect to the diagonal 
part of the Hamiltonian, 



H c — H\ + H2 + H12 



(E3) 



with 
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and where the marginal interband interaction reads 
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IKl. 
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Ko s K, 
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From there, we derive the RG equation for B xy ^ in a real 
space RG scheme^ 7 Parameterizing the running short 
distance cutoff as a(b) = ba, we obtain the flow of B xyi 
as 



dB. 



d\og(b) 



(1 - 9i) B 



xy,i • 



(E7) 
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Appendix F: Self-consistent derivation of the RKKY 
interaction 

In Appendix [El we established that the intermediate 
fixed point c of Fig. reached after integrating the flow 
of the relevant Overhauser fields, can be described as 2 
gapped and two gapless spinless Luttinger liquids. The 
remaining ungapped modes now mediate an RKKY 
interaction similar to case discussed in Appendix [C] We 
derive the latter in a mean field approximation where the 
<f>i + fields are locked to the minima of the cosine-terms. 
As one could expect, Gaussian fluctuations around this 
mean field have been shown to be unimportant^ and 
will be neglected here. We also neglect renormalizations 
of the Luttinger liquid parameters and effective velocities 
during the flow between the fixed points a and c of Fig. [7] 
since the gaps B xy ^ are strongly RG relevant and the 
associated flow is thus short. In addition, interband in- 
teractions are not important unless vp\ ~ vf 2 , see main 
text. In this approximation, the mean field Hamiltonian 
of the remaining gapless <pi- modes reads 



i=l,2 ' 

dz 4J7 e flE 



A 1 1 V 

{dAx-) , (Fl) 



2ir tt 

where the effective interaction between the 4*i- -modes is 
given by 



U eS = U 



Kic K 2c 



K Xs K x K 2s K 2 ■ 



(F2) 



This Hamiltonian can be diagonalized with a transfor- 
mation similar to Eq. (|B4I) . namely 
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(F4) 



where we used Su 2 2 - 



u 2 _ and with 



A 



V ( <5u i2-) 2 + {4U oB /tt) 2 U1-U2- + Su 2 12 _ 



(F5) 



This transformation brings the gapless sector of the 
Hamiltonian in the diagonal form 



k=± n J 



(F6) 



which in turn allows us to evaluate the residual RKKY 
interactions due to the gapless electronic modes 0+ and 
</>_ . As before, the residual RKKY interactions are deter- 
mined by the residual spin susceptibility. We start from 
Eq. (|C3j) . apply the transformation (|E2[) and drop the 
gapped fields <j)i + (which are simple constants in imagi- 
nary time and space and thus drop out). For the x and 
y components of the residual spin susceptibilities in the 
two bands, we obtain 
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4(27ra) 
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Evaluating these expressions and neglecting the velocity 
difference between <f> + and 0_ yields the spin susceptibil- 
ity as 



L sin (irg'i J / /3u.^\ 2 2s " 



2na 



(F10) 



17 



where T is the standard Gamma function and with 



which transforms the RKKY exchange coupling to 



2<4 



B 



(Fll) 



An analogous calculation can also be performed to calcu- 
late g[ 2 - Due to the uncertainty principle, an ordering of 
regularizes terms that depend on 9i + 1 13 i 28 i 30 Roughly 



•->i+ 



speaking, an expectation value involving the fields 6i + 
averages to zero. For the same reason, also the z compo- 
nents of the susceptibilities are non-divergent. 



Appendix G: Tyablikov decoupling and ordering 
temperatures 

Neglecting the magnetic field due to the electron po- 
larization, the symmetric sector of the nuclear spins is 
governed by the Hamiltonian 



J s rot -'(z. t - z s ) 



A-ij B-ij 

Bji Aij 






J s zz '(z t - 



with 



(G6) 



+A ji = Jf xl (zi - Zj) cos(2/c F i(zi - zj)) , (G7a) 
-B^ = J* x '{ Zi - zj) sm(2k F1 (zi - zj)) . (G7b) 



In this basis, the Overhauser field resulting from the 
Fourier components close to 2k p\ reads 
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and 



{zi), we obtain the real space 
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This Hamiltonian allows us to calculate, for instance, the 
critical temperature T* for the ordering of the nuclear 
spins in a helix of momentum 2k F i in the symmetric 
channel. It is at first useful to realize that the nuclear 
spin order as well as the RKKY interaction essentially 
only have Fourier components either close to ±2fcp i or 
±2kp2- Since the different minima of the RKKY interac- 
tion do not overlap, these two sets of Fourier components 
can be analyzed independently. We can thus rewrite the 
Hamiltonian as 
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Next, we apply the so-called Tyablikov decoupling which 
is based on the decoupling of the equation of motion of 
the nuclear spin Green's function. The latter is defined 
as 



G-V (*, = ~m ( [tfi„ (*), (*')J ) (G9) 

with I - (t) = (t) ± ?7 s z ;i|| (t). The equation of motion 
for this Green's function reads 



id t Gg* + _ (*, *0 = *(« ( £i„ (*). Ci„ (0 ) (G10) 
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Fourier transformation of the equation of motion finally 
yields 



d(t - t') e Mt-f ) (i dt g*£ + _ ( t) (Gil) 



u>G% + _(u) + (I+ iu (t),I-. i (t) ) 



(( H' RKK YiI s . 



S,»|| V " S,»|| 

re t. 



,1, 



*:)» 



with 



and treat the Fourier modes close to ikpi and ±fc_F2 
separately. Retaining at first only momenta q ~ ±2kpi, 
we can go to a rotated frame of reference where this or- 
der corresponds to a ferromagnetic polarization. This is 
achieved by the transformation 



cos{2kpxZi) sin(2fci?iZi) 0^ 
— sin(2fc^iZj) cos(2fcpiZi) 
1, 
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The commutation of the nuclear spin raising operator 
and the Hamiltonian can be evaluated as 
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The essential approximation of the Tyablikov decoupling 
now consists of the following simplifications, 
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ret. 

J i 

(G14b) 
(G14c) 



where [(•), (•)] + is the anticommutator, while [(•), (•)] de- 
notes the anticommutator as before. At this point, these 
decouplings are mathematically not justified. They have 
however successfully been applied before and in this sense 
benefit from an a posteriori justification ! 31 ' 32 Physically, 
the Tyablikov decouplings are indeed plausible. In the 
ordered state, the fluctuations around the ground state 
are small, such that one may approximate Z^v — > 
= const., which implies the first and second line of 
Eq. (|G14|) . The third line may be interpreted as fol- 
lowing from the fact that any decoupling would involve 
the ground state expectation value of a spin raising or 
lowering operator, which should, to good approximation, 
vanish in the semiclassical ground state. In addition, we 
also drop the purely local terms ~ An, J* (0). The equa- 
tion of motion for the nuclear spin Green's function is 
then given by 
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with Q = u> + i + . Now performing a second Fourier 
transformation to momentum space, defined as 



G^_(Cj) = 1 e iq{Zi - Zi) Gij,+-(") , (G16) 
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ret. 



The Green's function G™*-" (£>), i.e. the Fourier trans- 

form of G$L_ (t, t>) = -i6(t-t') ( [i^ (t), 1^ (f )] >, can 
similarly be shown to obey the equation of motion 

wG rct :(g ) w)= (G18) 



2A(q = 0)-A(q) -Ji(q) ^ 



Nl 



Ml) - J S Z (0) r*« 
N 2 , 



GT-(q^) 



We may thus define 



/TXX 2A(q = 0)-A(q)-J*(q) 
a q = (l s ) ^ > (G19a) 



Nl 



= (n A(q) - J s z (q) 



Nl 



(G19b) 



and obtain by plugging Eq. (|G17I) into Eq. (|G18I) that 



2 m 



2(7f) (a q -yJa*-Pl) 

& - (&l - /3 2 ) 

(G20) 

In a final step, we now have to self-consistently determine 
the magnetization (If). To this end, we realize that quite 
generally the Green's function and the associated density 
fulfill 

= J dun B {u) (-~Im{Gf_:( g ,w)}) , 

(G21) 

where ns is the Bose distribution function. From there, 
we find by Fourier transformation to real space that 



\ S,l 8,1/ J\J- / j 



Q2 v 

q 

(G22) 



« 2 - ?l 
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At the same time, we know from the angular momentum 
algebra that = (IAi)(/JV x +l)-I£-(i£) 2 . This 
leads us to a self-consistent equation for (I*) (through a q 
and fiq), 



pIN ± (IN ± + l) 



1 + 



■■n B {\/a 2 



P 2 q ) 



(G23) 



zation around the ordered state. It is thus given by 
(to) 2 ) ~ IN±(IN± + l)/3, since the system is dis- 
ordered for T = T^, and spin rotation symmetry is not 
broken at high temperatures. Because we are considering 
a finite size wire, the average ((-ff j) ) is also not allowed 
to jump at the transition. This sets p « 2/3. With this 
approximation, we can finally expand the Bose function 
for T = T c \ — ST in the small magnetization (If) and 
solve the self-consistent equation. Approximating finally 
IN±(IN ± + !)->■ {IN ± ) 2 yields 



where P IN ± (IN X + 1) = IN± (IN± + 1) - ((/^) ). 

Due to the appearance of ((/f^) 2 ), this equation is not 
sufficient to determine the magnetization for IN± > 1/2. 
One rather has to construct a set of 21 N± — 1 cou- 
pled equations by evaluating Green's functions of the 
form G" - ([I+,(I*) n I-]), see Ref. [p. To deter- 
mine the critical temperature, however, one may recall 
that for T = T c i, the magnetization (I*) considered 
here (i.e. the Fourier components close to q = 2fc.Fi) 
vanishes. For T ps T c \, the expectation value ((Zfj) ) 
is dominated by the large fluctuations of the magneti- 



2I 2 fa f . a'{q) 



with 



dq — 



3 \tt J &{q) 2 -P'{qf 



(G24) 



a'(q) = [2A(q = 0) - .4(g) - J*(q)] , (G25a) 
P\q) = [A(q) J*(q)] . (G25b) 

The critical temperature of the second band can be 
derived analogously. 
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